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ABSTRACT 



O 



We present a rigorous, detailed study of the generic, quantitative properties 
of gravitational lensing near cusp catastrophes. Concentrating on the case when 



O ! the individual images are unresolved, we derive explicit formulas for the total 

magnification and centroid of the images created for sources outside, on, and 
inside the cusped caustic. We obtain new results on how the image magnifications 
scale with respect to separation from the cusped caustic for arbitrary source 
CN \ positions. Along the axis of symmetry of the cusp, the total magnification \x 

^ ! scales as \i oc where u is the distance of the source from the cusp, whereas 

perpendicular to this axis, /i oc u~ 2 ^ 3 . When the source passes through a point 
\Q [ tto on a fold arc abutting the cusp, the image centroid has a jump discontinuity; 

we present a formula for the size of the jump in terms of the local derivatives of 
the lens potential and show that the magnitude of the jump scales as | ^5 1 1//2 for 
\u1\ <C 1, where is the horizontal distance between u and the cusp. The total 
magnifications for a small extended source located both on, and perpendicular to, 
the axis of symmetry are also derived, for both uniform and limb darkened surface 
brightness profiles. We find that the difference in magnification between a finite 
^ \ and point source is < 5% for separations of > 2.5 source radii from the cusp point, 

while the effect of limb-darkening is < 1% in the same range. Our predictions 
. ^ | for the astrometric and photometric behavior of both pointlike and finite sources 

passing near a cusp are illustrated and verified using numerical simulations of 
the cusp-crossing Galactic binary-lens event MACHO-1997-BUL-28. Our results 
can be applied to any microlensing system with cusp caustics, including Galactic 
binary lenses and quasar microlensing; we discuss several possible applications of 
our results to these topics. 

Subject headings: astrometry — stars: binaries, fundamental parameters — gravitational 
lensing 
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1. 



Introduction 



Over the past twenty years, gravitational lensing has grown from a mere curiosity to an 
important component of a large and diverse set of fields in astronomy. Its ubiquity is due at 
least in part to the fact that its effects are observable over a wide range of scales. This has 
enabled astronomers to use lensing to study everything from the smallest compact objects, to 
the largest structures in the universe, and almost everything in between. Despite the diversity 
of applications of gravitational lensing, the mathematical description of the phenomenon 
itself is both relatively tractable, and universal. In almost all cases, gravitational lensing 
can be described by a two-dimensional mapping from a lens plane to light source plane. 
Once this mapping is specified, all of the properties of a gravitational lens can be derived in 
principle. The observable properties of lensing, however, depend on the phenomenon to which 
it is applied. Therefore, lensing is traditionally divided into an number of different regimes, 
which are delimited by the observables. For example, the term microlensing is typically 
applied to the case when multiple images occur, but are not resolved. When multiple images 
are created by a gravitational lens, the separation between these images is typically of order 
the Einstein ring radius, 



where M is the mass of the lens, D = D 0S D i/ Di s , and D os , D a i, and Di s are the distances 
from observer to source, observer to lens, and lens to source, respectively. Thus, the term 
microlensing is applied when 6 E is less than the resolution. In this case, all one can measure 
is the collective behavior of all the images created by the lens, i.e., the total magnification, 
and the position of the center-of-light (centroid) of the image. In fact, if the observer, source, 
and lens were not in relative motion, than the individual image magnifications and positions 
would be fixed, and even these properties would not be measurable. The relative positions 
of the observer, lens, and source, and thus the magnification and centroid, are expected to 
vary on time scales of order the Einstein ring crossing time, 



where v± is the transverse speed of the lens relative to the observer-source line-of-sight. 
Fortunately, in the two regimes where microlensing has been discussed, the typical values of 
6 E , D i, and v± result in reasonable time scales. Typical values for the lens mass M, relative 
source-lens distance D, transverse velocity v±, and the resulting typical values for 9 E and 
t E for both Local Group and cosmological microlensing are given in Table 1. Also shown 
are typical values for the radius of the source's emission region, this radius in units of 
# E , p* = 0*/0e, and the time it takes the lens to cross the source, = p*t E . These latter 
parameters will be relevant to the discussion of finite source effects in §2.4. For the Local 
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M 


D 




v± 


t E 




P* 


t * 


Local Group 


1M 


10 kpc 


lmas 


100 km s" 1 


100 days 


lyuas 


10~ 3 


0.1 days 


Cosmological 


1M 


1 Gpc 


3 pas 


500 km s" 1 


10 years 


0.1/ias 


0.03 


100 days 



Table 1 Typical Microlensing Parameters. 



Group, a typical Einstein radius crossing time is t E = 0(100 days), whereas for cosmological 
microlensing, t E = 0(10 years). 

Of exceptional importance in microlensing is the existence of caustics: positions in the 
light source plane corresponding to the critical values of the lens mapping. On caustics, 
at least one image is formed that has formally infinite magnification (for a point-source). 
When a source crosses a caustic, both the total magnification and centroid of all the im- 
ages exhibit instantaneous, discontinuous jumps. These jumps are averaged out over the 
finite source size; however, it is generically true that large gradients in the magnification 
and centroid exist near caustics. Furthermore, microlensing caustics have several important 
and useful properties. First, the large magnification results in a large photon flux from the 
source. Second, the large gradient in the magnification and centroid with respect to source 
position effectively implies high angular resolution. Finally, the highly localized nature of 
the high-magnification and large centroid-shift regions created by caustics results in charac- 
teristic, and easily-recognizable features, in both astrometric and photometric microlensing 
curves. Many authors have suggested exploiting these properties of caustics to study a num- 
ber of astrophysical applications, i.e., stellar multiplicity (Mao & Paczyhski 1991), stellar 
atmospheres (Gould 2001), individual microlens mass measurements (Graff & Gould 2002), 
microlens mass functions (Wyithe, Webster, & Turner 2000), properties of the emission re- 
gions of quasars (Wambsganss, Paczynski, & Schneider 1990; Agol & Krolik 1999; Fluke & 
Webster 1999; Wyithe, Webster, & Turner 2000), and lens transverse velocities (Wyithe, 
Webster, & Turner 1999). See Gaudi & Petters (2002; hereafter Paper I) for a more thor- 
ough discussion of the uses of astrometric and photometric microlensing observations in the 
presence of caustics in both Local Group and cosmological contexts. 

Although the caustic curves of gravitational lenses exhibit an enormously rich and di- 
verse range of properties, it can be shown rigorously that each stable lensing map has only 
two types of caustic singularities: folds and cusps (Petters, Levine, & Wambsganss 2001, p. 
294). Each of these two types of singularities have generic and universal properties and, in 
particular, each can be described by a polynomial mapping from the lens plane to light source 
plane. The coefficients of these mappings depend on local derivatives of the dimensionless 
surface potential of the lens. In Paper I, we used the mapping for a fold singularity to derive 
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the observable properties of gravitational lensing near folds, paying particular attention to 
the case of microlensing, in which the images are unresolved. We derived analytic expressions 
for the total magnification and centroid shift near a generic, parabolic, fold caustic. We then 
showed how these expressions reduce to those for the more familiar linear fold, which lenses a 
nearby source into two equal magnification, opposite parity images whose total magnification 
is proportional to u^ 1 ^ 2 , where u is the distance of the source to the fold caustic. We then 
generalized these results to finite source sizes. Finally, we compared our analytic results 
to numerical simulations of the Galactic binary-lens event OGLE-1999-BUL-23, in which 
the source was observed to cross a fold caustic. We found excellent qualitative agreement 
between our analytic and semi-analytic expressions for the photometric and astrometric be- 
havior near a fold caustic, and our detailed numerical simulations of the second fold caustic 
crossing of OGLE-1998-BUL-23. 

In this paper, we present a similarly detailed study of the generic, quantitative properties 
of microlensing near cusps. Although fold caustic crossings are expected and observed to 
dominate the sample of caustic crossings in Galactic binary events (Gaudi & Gould 1999; 
Gaudi, Graff, & Han 2002; Alcock et al. 2000), cusp crossings will nevertheless represent 
a non- negligible fraction of all caustic crossing events. In fact, at least two cusp crossing 
events have already been observed, the Galactic bulge events MACHO 97-BLG-28 (Albrow 
et al. 1999a), and MACHO-1997-BLG-41 (Alcock et al. 2000; Albrow et al. 2000). It is 
interesting to note that the analyses and modeling of these events were performed using 
entirely numerical methods. As we discuss in some detail (see §5), we believe that the 
analytic results derived here are particularly amenable to the analysis of MACHO 97-BLG- 
28 and similar events. In the cosmological context, the role of cusps versus folds is less clear, 
due primarily to the more complicated structure of the caustics themselves. However, it has 
been shown that, in the limit of high magnifications, the inclusion of cusps alters the form 
of the total microimage magnification cross section, due to the lobe of high-magnification 
close to and outside the cusped caustic (Schneider & Weiss 1992). Lensing near cusps has 
been substantially less well-studied than lensing near folds, and as a result, useful analytic 
expressions for the observable properties are few. Previous studies have focused almost 
exclusively on the magnification of the images created by a cusp singularity (Schneider & 
Weiss 1992; Mao 1992; Zakharov 1995, 1999). Here we study all the observable properties 
of gravitational microlensing near cusps, including the photometric (total magnification) 
and astrometric (light centroid) behavior, for both point sources and extended sources with 
arbitrary surface brightness profiles. 

The layout of this paper is as follows. In §2, we derive analytic expressions for the image 
positions, magnification, and light centroid for sources near a generic cusp. In §2.1 we define 
the observable microlensing properties. In §2.2, we start with the generic expression for the 
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mapping near a cusp, and derive all the properties for the local images for sources exterior 
to (§2.2.2), on (§2.2.3), and interior to (§2.2.4) the caustic. We generalize the discussion to 
include images not associated with the cusp in §2.3, and study extended sources in §2.4. In 
§3, we illustrate the observable behavior near a cusp by numerically simulating the Galactic 
binary-lens cusp-crossing event MACHO 97-BLG-28, and directly compare these numerical 
results with our analytic expressions in §4. In §5, we suggest several possible applications 
in both local group and cosmological contexts. We summarize and conclude in §6. We 
note that, for the sake of completeness, we include some results that have been presented 
elsewhere. Combined with the results from Paper I, the results presented here describe the 
observable properties of gravitational microlensing near all stable singularities. 



2. Analytical Treatment 

2.1. General Equations for Magnification and Image Centroid 

The lens equation due to a gravitational lens with potential tp is given in dimensionless 
form as follows: 

u = 0-a(0). (3) 

Here 6 = r/(0-ED o i) and u = s/(0eD os ), where r and s denote the proper vector positions 
in the lens and light source planes, respectively, and ct(0) = V^(0). Note that if k is the 
surface mass density of the lens in units of the critical density S cr = (c 2 D 0S )/(4nGD iDi s ), 
then the gravitational lens potential ip is related to k by 

rj}{e) = - f d0'K{0')\n\0-0'\. (4) 
* Jr 2 

The lens equation induces a lensing map rj from the lens plane L into the light source plane 
S defined by rj(6) = — ot(0). The critical curves of t] are defined to be the locus of all 
critical points of rj, i.e., the set of solutions C of the equation 

det[Arj](e c ) = 0, (5) 

where Arj is the Jacobian matrix of 77. The caustics of r\ are its critical values, i.e., the set 
of all points u c = r](0 c ). 

For a fixed source position u, the solutions in L of the lens equation determine the 
lensed images of u. Suppose that a source at u has N{u) images. The magnification of the 
ith image 0; = 0i(u) is given by 
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Then the total magnification is 



N(U) 



(7) 



i=l 



For a source at u, the image centroid or center- of -light, denoted by c i(w), and shifted image 
centroid, which is denoted by 56 c \, are defined respectively as follows: 



The image centroid is defined analogously to a center of mass. However, the image centroid 
is much more complicated since images can suddenly be created or annihilated as the source 
moves over caustics. For the case of gravitational microlensing, the total magnification and 
shifted centroid are observables, with the latter expected to be accessible to, e.g. the Space 
Inter ferometry Mission (SIM). 



A rigorous analytical study of the total magnification (photometry) and image centroids 
(astrometry) for most physically reasonable lens models (e.g., binary point mass lenses) 
is very difficult. Our goal here, as in Paper I, is to study the local, generic behavior of 
these observables near an arbitrary, but fixed, critical point of a given type. In Paper I, we 
studied the behavior near a fold singularity; here we study the behavior near a cusp. In 
order to study the local behavior near singularities, we Taylor expand the scalar potential 
tp in the neighborhood of the singularity, keeping terms to a specified order in the image 
position 0. We note that the order at which one chooses to truncate the expansion is 
somewhat arbitrary; higher accuracy can always be achieved by keeping terms of higher 
order (at the expense of higher complexity, of course). However, it is generically true that 
fold critical points are characterized by partial derivatives of the lensing map up to second 
order, while cusps are classified by partials up to third order (Petters et al. 2001, pp. 370- 
371). As a result, the Taylor expansion of the potential is usually taken up to second order 
for the fold and third order for the cusp. See Schneider et al. (1992, pp. 187-188) and 
particularly Petters et al. (2001, pp. 341-353) for more discussion on truncation of Taylor 
expansions of the potential. After introducing an orthogonal change in coordinates, one can 
thus obtain a simpler equation that approximates the quantitative, generic behavior of the 
total magnification and image centroid of sources near a critical point. Consider a generic 
lensing map r] whose caustics have cusps. In such a case, there are at least two cusps since 




N(U) 




(8) 



2.2. Local Case: Point Sources Near a Cusp 
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the total number of cusps is always even (see Petters 1995, p. 4282, for a cusp-counting 
formula). In this section we will fix one of the cusps and study the behavior of rj in the 
vicinity of that cusp. 



2.2.1. Local Case: Magnification and Centroid of Images Associated with the Cusp 

Without loss of generality, we can assume that the coordinates 6 on the lens plane L 
and u on the light source plane S are translated so that the cusp lies at the origin of S 
and 77(0) = 0. There is an orthogonal change of the coordinates 6 and u which is the same 
in the lens and light source planes such that the lensing map r\ can be approximated by the 
following simpler mapping in a neighborhood of the origin (Petters et al. 2001, pp. 341-353; 
Schneider et al. 1992, p. 193): 

U\ = c6 1 + h -Q\, u 2 = 60102 + a6l (9) 

where 

a = -^2222(0), b = -^122(0)^0, c = l-^ii(0)^0, 2ac-&V0. (10) 

Note that in (9) we are abusing notation somewhat by using the notation = {61,62) and 
u = (1x1,1x2) for the global lens equation (3) to express the local lens equation. For the 
explicit relationship between the local coordinates in equation (9) and the global ones in 
equation (3), see pp. 344-345 of Petters et al. (2001). It is also important to add that 
the partials in (10) are with respect to the original global coordinates of the lens equation. 
Figure 1 illustrates the coordinate systems and basic properties of lensing near a cusp. For 
this figure, we have adopted the coefficients a = 1.69, b = 2.45, and c = 2.0. These are the 
appropriate values for the cusp in the observed event MACHO 97-BLG-28. See §4. 

We shall now use the lens mapping determined by equation (9) to study the approx- 
imate behavior of the magnification and image centroid near a cusped caustic curve. The 
magnification matrix A for the cusp is the Jacobian matrix of equation (9), 

A= (c bd 2 

V b9 2 &0i + 3a9 2 2 

which has the determinant 

det A = bc9 1 + (3ac - b 2 )Q\. (12) 

The cusp critical point of the mapping in equation (9) is at the origin of the lens plane and 
mapped to a cusp caustic point at the origin of the light source plane. By equation (12), the 
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critical curve is a parabola, 

ft = (13) 
Equations (9) and (13) yield that points u c = (i£ c ,i, 1^,2) on the caustic are given by 

3(b 2 - 2ac) n2 b 2 -2ac n3 

u c ,i = -^-^ — '-el u cfl = — - — e 3 2 , (14) 

where 6 2 acts as a parameter along the caustic. The magnitude of the curvature of the 
caustic is 

" L cA "c,2 — "c,2 "c,ll 



k 



K.1 W'n-U'n U" A \ \b C\ 1 



[KJ 2 + « 2 ) 2 ] 3 / 2 3(b 2 + c 2 ) \b 2 - 2ac| \9 2 \ ' 

where the primes indicate differentiation with respect to 9 2 . The curvature diverges as the 
cusp is approached (9 2 — > 0). 

By equation (14), the caustic can also be described as the level curve 

8b 3 

C( "^ 27c*(2 a c-^ + "' = - (15) 

The origin is a positive cusp if 2ac > b 2 and a negative cusp if 2ac < b 2 . Without loss of 
generality, we shall assume that the cusp is positive. For this case, both a and c have the 
same sign. The conditions C(u) > 0, C(u) = 0, and C(u) < determine points u outside, 
on, and inside the cusped caustic curve, respectively. 

Let us determine the local lensed images and their magnifications for sources near the 
cusped caustic. The expression for u\ in equation (9) yields 

* = *-±* (i.) 

Plugging equation (16) into the expression for u 2 in equation (9), we obtain a cubic equation 
for 6 2 in terms of the source position {u\,U2)\ 

/(£ 2 )^ 2 3 + p# 2 + q = 0, (17) 

where 

P = PW ^ U1 ' q = qM = -2^ M2 - (18) 
In other words, the lensed images are points 6 = (61,62) with 6\ given by equation (16) and 
6 2 a real root of equation (17). The magnification of each lensed image of u is then 

"< e > = raor (19) 
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where 

det A = bu\ + — (2ac — b 2 )9 2 (20) 

with 2 a real zero of /. 

For a source at angular vector position u, let 0*{u) = {6\{u) , O^iu)) be any lensed 
image of u associated with the cusp, i.e., 

*.(«)= (^-^[^2»] 2 ,^(«)), (21) 

where ^(u) is a solution of cubic equation (17). Then equation (14) yields that 

_ 3(b 2 -2ac) 2 _ fe 2 -2ac 3 
M ci 26 ^ 2 ' ' c ' 2 ~~ c ^ 2 ' ' ^ ' 

is a point on the caustic. It is important to note that even for a source position u = u c at 
a caustic point, there is no guarantee that the caustic point u* c = (u* i,u* 2 ) equals u c (see 
Eq. [53] below). Now, consider the tangent line to the cusp, which happens to coincide with 
the Mi-axis (due to our coordinate choice). We shall also refer to this tangent line as the 
axis of the cusp, since it coincides with the axis of symmetry of the cusp. Note that by (11) 
the magnification matrix A(0) at the origin, maps the entire lens plane into the previous 
tangent line. Define the horizontal distance between u = (^1,^2) and u* = {u* cl ,u* c2 ) to be 
the distance between u and u* along the direction parallel to the tangent line to the cusp. 
Denoting the horizontal distance by Au^ h , we obtain 

A«; >h = k - <,| = k - 3(62 ~ 2flc) (6^) 2 |. (23) 

By equation (20), we see that the magnification of 0* is inversely proportional to the hori- 
zontal distance between u and u* : 

^ = (24) 

We shall see that if the source is on the Mi-axis, then the magnification of the image 0* (only 
one image occurs locally for sources outside the caustic) is proportional to u^ 1 , where u is the 
distance between the source and cusp (see Eqs.[36] and [70] below). If the source is on the 
w 2 -axis, then the magnification of 0* (again, only one image occurs locally) is proportional 
to u~ 2 l z (see Eqs. 33). In both cases, the horizontal distance coincides with the ordinary 
distance. 



We note that in Paper I, we showed that the magnification of an image 6* associated 
with the fold caustic is inversely proportional to the square root of the perpendicular distance 
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of the source from the fold. For a generic local fold lensing map that maps a fold critical 
point at the origin to a fold caustic point at the origin, the perpendicular distance is defined 
relative to the direction perpendicular to the tangent line to the fold caustic at the origin. 
Note that the magnification matrix A(0) for the fold also maps the entire lens plane into the 
tangent line. 

Now, the reals zeros of the polynomial / in equation (17) are characterized using the 
discriminant of /, i.e., 

B W-(I) 3+ (I) 2 =(2^ C W- < 25 > 

By equation (25), the conditions D{u) > 0, D{u) = 0, and D(u) < also determine, 
respectively, points u outside, on, and inside the cusped caustic curve. It will be seen that 
for a source outside, on, or inside the cusped caustic, there are locally one, two, or three 
images, respectively. These images lie in the neighborhood of the origin, which is a point on 
the critical curve. The magnification and centroid of the images associated with the cusp 
will be denoted by /ii oc and 0\ oc , respectively, and called the local magnification and local 
image centroid. We shall now determine fi\ oc and 0\ oc for source positions near the cusp that 
are outside, on, and inside the caustic — see equations (80) and (81) for a summary. 



2.2.2. Local Case: Source Outside the Cusp 

For a source position u outside and near the cusp (i.e., D(u) > 0), equation (17) has 
one real root, namely, 



0?\u) = (-| + v ^)) 5 +(-|- v ^))% 

u 2 + 



2ac-b 2 

Consequently, there is one lensed image (locally) 

JOUt 



o (»= (^-^rNi 2 , r («) ! 



(26) 
(27) 

(28) 



To determine the parity of OU (u), it suffices to find the sign of det A(0 out (u)) for a single 
source position outside the caustic. A source lying on the w 2 -axis either above or below 
the origin is outside the caustic curve (since (7(0,1^) = u\ > for w 2 ^ 0). Furthermore, 
#° ut (0,w 2 ) = (-q)i, which is nonzero for w 2 7^ 0. It follows that 



detA[0 out (O, M2 )] = -(2ac-6 2 )(-q)l > 0, 



(29) 
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where we used the fact that the cusp is positive (i.e., 2ac — b 2 > 0). Hence, the image 9 ont 
has positive parity. In the case of lensing by n-point masses, there are no maximum images 
(Petters 1992) and so out is a minimum image. 

The magnification of the image out (u) is 



K o u c >)=M0 out N) = 



\b\Au° u ^ 



where 



Au°* = \u 1 --(b 2 -2ac)[er(u) 



(30) 



(31) 



As u approaches the cusp at from outside the caustic (i.e., the constraint D{u) > 
holds for m ^ 0), we have 0§ ut (it) — > 0. Consequently, the magnification fi°^ continuously 
increases without bound as the cusp is approached from outside. This is unlike the case 
for a source crosses a fold caustic since then a sudden, infinite, discontinuous- jump occurs 
in the magnification (for point sources). The continuous increase in is due to a lobe 
of high-magnification outside the cusp (see Fig. 4, and also pp. 334-335, 480-485 of Petters 
et al. 2001). 



We now consider the magnification of sources on the axes, exterior to the cusp. For 
sources on the -u 2 -axis and outside the cusp, we have p = and (thus) C(u) = u 
resulting image position is therefore 



% The 



out (0,u 2 ) = 



2c 



2c \ 2ac - - '"J 



\ 2 /3 



U-2 



2c 



2ac - b 2 



u 2 



1/3 



(32) 



and the magnification is given by, 

Ko u c(0,« 2 ) =M0° Ut (0,« 2 )) 



Ur2 
U 2 



2/3 



(33) 



where we have defined the characteristic rise scale u r2 for trajectories along the w 2 -axis, 



u r2 = 



27 



c 2 (2ac - b 2 ) 



-1/2 



(34) 



Thus, the magnification diverges as u 2 2 ^ 3 for u 2 — > 0. 



If the source is on the Mi-axis and outside the cusp (i.e., U\ > 0), then q = 0. Equation 
(26) yields #2 Ut (wi,0) — 0, so the image position is simply, 



6 



out 



(«!,0)= (^,0), 



(35) 
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and the magnification is given by, 



/C(o,«2)=M0 out K,o)) 



(36) 



where we have defined the rise scale u r ±, 




(37) 



Hence, the magnification diverges as u t for u\ — > 0. 

For a source moving along a path u(t) (not necessarily rectilinear) lying outside the 
cusp, the local image centroid follows the motion of the lensed image OU (u(t)) (see the 
path of out in Figure 1), i.e., 



Observe that 0°^(0) = 0. Figures 1 (b) and (f) illustrate the trajectories of the image and 
centroid, respectively, for a source in rectilinear motion. 

Let us now consider the behavior of the lensed image 6 out when the source at u ap- 
proaches the caustic. First, as u approaches the cusp from outside, we obtain 6 out (u) — > 0. 
Second, at a point Uq = (it®, m°) on one of the fold arcs meeting the cusp, equation (26) shows 
that #2 Ut (™o) = 2(-q /2) 1/3 . Since D(u ) = 0, it follows from equation (25) that 



where p = p(w?) and q = q^). Note that equation (39) implies — p > 0, which yields 



(since 2ac — b 2 > 0). Hence, the limiting position of 8 out (u) as u — > u from outside the 
caustic is 



(38) 




(39) 



bu\ < 



(40) 




(41) 



2.2.3. Local Case: Source on the Caustic 



Assume that the source is on the caustic (i.e., D{u) 
point u = 0, then 

p = q = 0, 



= 0). If the source is at the cusp 

(42) 
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and equation (17) has a real triple root, 2 (O) = 0. This yields an infinitely magnified lensed 
image at the origin of the lens plane. Consequently, for a source on the cusp, the local image 
centroid is at the origin (as can be seen from Eq. [38]): 

= out (O) = 0° o n c (O). (43) 



If the source is at a point u = (u^u®) on one of the fold arcs abutting the cusp (i.e., 
D{uq) = and u ^ 0), then 

Poqo ^ 0, (44) 
so equation (17) has two real roots, one a simple zero and the other a double zero: 

er-M = 2 (=£ V , er ' d M = - J . (45) 



The corresponding lensed images are as follows (since D(u ) = 0): 

c 



6 




f-<io 


6c 


Po,- 


V 2 


2b 


Po,2( 




3c 


, 2 , 



(46) 



fold '>o) - + g Po, 2 (^J J = 0° ut (n o ). (47) 

Note that fold ' d (O) = = fold,s (O) and the rightmost equality in equation (47) is a conse- 
quence of equation (41). 

Using equation (20), equation (45), and D(u ) = 0, we get 

det A{0 Md ' s {u )) = bu\ - 2{2ac - 6 2 )p = -36w° > 0, (48) 
detA(0 fold 'Vo)) = K--(2ac-6 2 )p = 0, (49) 

where the positive condition follows from equation (40). In other words, the image fold,d (u ) 
has infinite magnification and so is located on the critical curve. On the other hand, the 
image 6 iold ' s (u ) is finitely magnified and has positive parity: 

^ loc(e -.. (Uo)) = _L_ = _±, (50) 

where the rightmost equality follows from equation (40). Now, the root 6*2° ld ' s (wo) determines 
the caustic point *4° ld > s = (u™ d ' s , u^' s ) , where 



foid.s 3(6 2 — 2ac) m foid,s/ xi2 foid.s b 2 — 2ac iold s 3 
«ci = 2b ' 2 ^ °^ ' Uca = : ^ 2 ^ °^ ' ^ 
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By equation (45), we have 



fold,s 



4u° v (52) 



"c,l 

Consequently, the horizontal distance from the source at the fold point u to the caustic 
point u { ° ld,s is nonzero: 

Au[ ^ s = |u° - u { c f> s \ = 3|u?| ^ 0. (53) 
By equation (24) or (48), the magnification of 6 iold,s (u ) can be expressed as 

/4ot d ' S K) = (54) 
|6|A« ljH ' 

In other words, the magnification of fold ' s (u o ) is not inversely proportional to the horizontal 
(or ordinary) distance from u to the caustic (which is zero since the source sits on the caus- 
tic), rather from u to the caustic point n^ old ' s . Analogous to 2 old ' s (M o ), the root 6*2° ld ' d (uo) 
in equation (45) yields a caustic point i4° ld ' d = (u°^' d ,u°^' d ). Since 

fold.d 3(6 2 — 2aC) m fold,d/ \12 /rr\ 

u d = 2fo ^ 2 ^ =Ml ' ^ ■> 

the horizontal distance from u to u { ° ld ' d is zero, 

Au^ = |u° - <f d | = 0. (56) 

Since (J^ ld,d (u ) = l/(\b\ Au° ld,d ), the vanishing of the horizontal distance gives the earlier 
result that fold ' d (n,o) is infinitely magnified. 

As u approaches a fold point u from outside the caustic, we have D(u) — > and 
9% ut (u) — > ^2° lds (' u 'o) (using Eq. [27]). This implies that the outside local image centroid 
obeys 

<CH -> fold '> o ). (57) 

However, when u reaches u , the centroid 0°^( u o) n °t only meets 8 iold ' s (u ), but also the 
infinitely magnified image iold ' d (u ). The latter image suddenly appears when u = Uq and 
dominates fold ' s (u o ) 

in magnification. This causes the image centroid to have a double-value 
(see Eq.[79] below) and creates a discontinuous jump from iold ' s (u o ) to 6 iold ' d (u ). Further 
discussion of this jump discontinuity will be given below (see equation Section 2.3.2 and 
Figure 2). 
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2.2.4- Local Case: Source Inside the Cusp 



For sources inside and near the cusp (i.e., D(u) < 0), equation (17) has three real simple 
roots: 



C>) = 2 A /^cos 



-P 



6^» = -2^-^oos 



-p 



■d(u) — IT 



^» = -2^cos 



-P 



$(U) + 7T 



where 



•&(u) = cos 



-i 



(-q/2) 



(58) 



(59) 



.V / Fp73F. 

Since -D(u) < 0, it follows that (p/3) 3 < -(q/2) 2 < 0, which yields p < and bu x < (since 
2ac - b 2 > 0). In addition, we have \q/2\/y/(--p/3) 3 < 1, or, -1 < (-q/2)/ y/ ^(-p/3) 3 < 1. 
Consequently, 

< < 7T (60) 

for w inside the caustic. The values d(u) = and d(u) = n occur when u is on the top 
and bottom fold arcs, respectively, that abut the cusp (see discussion below). The value 
■&{u) = 7r/2 is for u on the Mi-axis inside the caustic. 

A source inside and near the cusp has three local images determined by equation (58): 

in » = f- + |W 

\ c 3c 

e^iu) = (^i + ^pcos 2 

c 3c 




(61) 



We saw that the condition D(u) < yields bui < and < < n for sources inside the 
caustic. The latter, implies 



detA(0' m '\u)) = &ui 
detA(O in ' 2 (u)) = 6ui 



1-4 cos 



1-4 cos 



>0, 



■d(u) — 7T 



>0, 
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detA(0' m ' 3 (u)) = 6ui 



I - I cos- ! ) < 0. 

(62) 

In other words, the lensed images and in ' 2 have positive parity, while 6 m ' 3 has negative 
parity. It is useful to note that 

det A{O in '\u)) = 0, for 6 = tt, 
detA(0 in ' 2 (u)) = 0, for = 0, 
detA(6» in ' 3 (u)) = 0, for# = 0,7r. 

(63) 

Let us now consider the behavior of three images in equation (61) near the caustic - 
see Figures 2 and 3. Equation (61) shows that as u — > (source approaches the cusp) from 
inside the caustic, we obtain 

6> in 'V)^0, £ = 1,2,3, (64) 

i.e., all three images merge at the point on the critical curve that maps to the cusp. We now 
investigate the image behavior as the source approaches a fold caustic. There are two fold 
arcs abutting the cusp, one in lying above the «i-axis and the other below. If the source is 
on one of these fold arcs, then D(u) = and pq ^ 0. At a point u on the fold arc above 
the Mi-axis, we have (— q /2) = a/ (— Po/3) 3 , which yields $(u ) = 0. By equation (61), as 
the source position u approaches u from the interior of the cusped caustic, we have 

d in '\u) — > fold ' s K), in ' 2 (u) — > fold ' d (u o ), in ' 3 (u) — ► fold ' d (u o ). (65) 

If the source crosses the fold transversely through the point Mo and continues on to outside 
the fold, then the lensed image moving along m ' 1 (u) passes through 6 iold ' s (u ) and continues 
along ont (u). In addition, the other two images moving along m,2 (w) and in,3 (u) merge 
into a single lensed image at fold ' d ('Uo), and then disappear as the source moves outside the 
caustic. These results are depicted in Figure 2(b). Similar results occur if the source travels 
from outside the caustic approaching a point Uq lying on the fold arc below the i^-axis, 
see Figure 3(b). In this case, (— q /2) = —a/ (— Po/3) 3 , which yields ~&(u ) = n. As the 
source approaches u Q from outside, there is one image locally and the image moves along 
out . Equation (61) shows that as u — > u from outside the caustic, we have 

6 out (u) — > fold ' s (u o ) = 0' m ' 2 (u ). (66) 

When the source reaches u , an infinitely magnified lensed image fold ' d ('Uo) suddenly ap- 
pears. As the source moves transversely through u and enters inside the caustic, equa- 
tion (61) yields that fold ' d (u o ) splits into two images that move along m,1 (u) and 6 m,3 (u). 
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Equivalently, if a source at u moves from inside the caustic to a point u on the fold arc 
below the tii-axis, we have that in>1 (u) and in,3 (ii) merge at fold ' d (tt o ): 

6 [n ' 2 (u) — > 6> fold ' s (u ), in ' 1 (u) — > fold ' d (uo), in ' 3 (u) — > fold ' d (u o ). (67) 



The magnifications of the images m ' (u) are 

fi(u) = M^W) = — * £ = 1, 2, 3, (68) 

|6|A Wl ^ 

where the horizontal distances are given by 

Aufj 1 = |«i (1 - 4cos 2 [t?(u)/3]) |, Aufj. 2 = |«i (1 - 4cos 2 [(t?(ii) - tt)/3]) |, 
Aufj. 3 = |ui (1 - 4 cos 2 [(0(u) + tt)/3]) |. 

(69) 

Observe that if the source is inside the caustic, but on the tii-axis, then i?(ui,0) = 7r/2, 
which yields 

0) = 2^ = 2^- ^ °) = 2^ = ^ ""(^ °> = w = sir- 

(70) 

In other words, the magnification is inversely proportional to the distance of the source (on 
the tii-axis) from the cusp. The magnification of the images are shown in Figures 2(b) and 
3(b) for a source in rectilinear motion. It follows from equation (70) that 

^(«i,0)=^(«i,0)+^(«i,0). (71) 

In fact, this relationship holds for any position u inside and near the cusp (e.g., Schneider 
& Weiss 1992; Zakharov 1995, 1999): 

/ 4»=/4»+/4», (72) 

where D(u) < 0. Hence, the total magnification of the three local images is 

,C( W ) = M» + ,4» + /4» = 2 [M» + /4»] . (73) 



The image centroid of the three images for source positions inside and near the cusp is 
given by 

^e^ 1 + /4 n in ' 2 + /4 n in ' 3 _ i 



^loc - 



^in(l,2) 0in,3 



(74) 
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where O^c' 2 ^ is the centroid of the images and m ' 2 , i.e., 

m(1 , 2) _ /^,i0*.i + 

^loc - 7^1 + • ^5) 

The image centroid for a source in rectilinear motion is shown in Figures 2(f) and 3(f). 

Equation (65) shows that as u — > uq from inside the caustic, where uq lies on the fold 
arc above the ^i-axis, we obtain // in,1 (it) — > /^"^(ito), where < {i[° ld ' s (u ) < an d 
fi m ' 2 (u) — > oo. This yields 

°> ..in.T^.N I .. in .2^ ^ 1 asw^uo- (76) 



//^(it) + yu in ' 2 (w) //^(w) + /i in ' 2 (n) 

Similarly, equation (67) shows that if u is on the fold arc below the lii-axis, then /t m ' 2 (it) — > 
/J^f ,s (u ) and fi m,1 (u) — > oo. Consequently, 

• as w -> u Q . (77) 



//^(it) + /j, in ' 2 (u) * n m > l (u) + /i in > 2 (u) 

It follows that as u — > it from inside the caustic, we have 0£< 1,2) (u) -> 0™ and in ' 3 (it) -> 
fold ' d . Hence, 

- e Md ' d ( Uo ). (78) 

as u — > u from inside the caustic. Equations (57) and (78) then show that for u = w 0) 
the local image centroid is double valued and so has a jump discontinuity (Figures 2(f) and 
3(f)): 

loc K) = {0 fold '> o ), fold 'Vo)}. (79) 



2.2.5. Local Case: Summary of Local Results 



The magnification of the local images associated with a source at angular vector position 
u near the cusp can now be summarized as follows: 



f /C(«) 



Vloc(u) = < OO, 



■u outside caustic, 



■u on caustic, 



(80) 



2 [//"(it) + //If (it)] w inside caustic, 



where the magnifications /t° u c * and /t^ n , £ = 1, 2, are given by equations (30) and (68), respec- 
tively. These magnifications are shown as a function of source position in Figures 1-4. The 
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image centroid is given by: 

' e out (u) 



01oc(«) = < 



) /jfoldjd/ \ /jfold,s/ \ I 
{^loc («0), loc («0)|, 

o, 



i/, outside caustic, 
i/, = uq on fold arc, 
it = on cusp, 



(1/2) 0£? ,2) (u) + in ' 3 (w)] u inside caustic, 



(81) 



where the angular vector ont is given by equation (28), #[°|, d ' d by equation (46), #[°|. d ' s 
by equation (47), 111 ' 3 by equation (61), and 0"^' 2 ^ by equation (75). Note that the above 
equations can be applied to a general (not necessarily rectilinear) source trajectory u = u(t). 
The local centroid is illustrated in Figures 1(f), 2(f), and 3(f) for rectilinear source motion. 



2.3. Global Case: Point Sources Near a Cusp 

2.3.1. Magnification and Centroid of All Images 
Suppose that the source moves rectilinearly, i.e., its trajectory is given by 

u(t) = u c + (t — t c )u, 



(82) 



where u c is the position at which the source intersects at time t = t c the tangent line to the 
cusp, which coincides with the the Mi-axis, and it is the constant angular velocity vector of 
the source, 

'cos0 sin s 



* = hr'-irj- (83) 

Recall that t& = D o i0e/v±, where 6>e is the angular Einstein radius and vj_ is the proper 
transverse speed of the lens relative to the observer-source line-of-sight. Here is the angle 
of the source's trajectory with respect to the Mi-axis (i.e., tangent line to the cusp). 

As the source moves along the rectilinear path in equation (82), the angular positions 
of the images unaffiliated with the cusp follow continuous trajectories and their magnifica- 
tions are positive, finite, and continuous. Moreover, we assume that for source trajectories 
near the cusp, the image magnifications and positions of the unaffiliated images are slowly 
varying functions of the source position u. Denote the centroid and total magnification of 
the unaffiliated images by #o and /iq, respectively. Then the magnification fi tot and image 
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centroid C \ of all the images are given by 

^tot(u) = A*ioc(w) + A*o(«) (84) 

and 

c \{u) = [//ioc(w)0i oc (w) + // o («)0o(u)] • (85) 

The photometric and astrometric observables /i to t and #i oc , respectively, can be expressed 
as functions of time by replacing u in equations (84) and (85) by the rectilinear source 
trajectory u(t) in equation (82), and Taylor expanding 6 and /i to first order about the 
crossing time t — t c of the tangent line to the cusp. This yields: 

IMot(t) = fJ,]oc(t) + Ho(t), (86) 

where H\ oc (t) = fii oc (u(t)) and fi (t) = /x , c + (t - £ c )/i(o,c) with /i (0 , c ) = dfi /dt\ t =t c , and 

Od(t) = hoc(t)6>loc(t) + ^o(t)6> (t)] , (87) 

Utot \t) 

where 0i oc (*) = O loc (u(t)) and O = ^o, c + (* - £ c )0(o, c ) with (O;C ) = d0 o /dt|t=t c . 



Image Centroid Jump-Discontinuity at Fold Crossing 

Assume that a source in rectilinear motion u(t) transversely crosses the point u on one 
of the fold arcs abutting the cusp. As u(t) — > w from inside (respectively, outside) the cusp, 
denote the limiting value of the image centroid O c \(u(t)) by 6^(u ) (respectively, 0^(u o )). 
The jump-discontinuity vector at u for the image centroid and shifted image centroid are 
the same: 

e nP K) = c >o) - *d(«o) = ^5 mP («o). (88) 
As u — > u from inside the cusp, we have 

P'n(u) 1 . , 

0, (89) 



/itot(w) 1 + A t ioc(wo)/A t o(wo) 

(since fi\ oc (u) — > oo and Ho{u) is generally non-divergent). Therefore, from equation (85), it 
follows that 0^(ito) = fold ' d ( w o); i- e - the contribution from the images not associated with 
the cusp is negligible. In contrast, as u — > u from outside the cusp, we have 

Mn) 1 ^0, (90) 



/i tot (w) 1 + /i° o u c t (u )//i (u ) 
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where the last inequality holds provided that u ^ 0. Thus, we generally have that 0^(i* o ) 7^ 
6 iold,s (u ), and the contribution from images not associated with the cusp can be substantial. 
However, due to the fact that the magnification of the image exterior to the cusp is ^ 1 
for positions close to the cusp, it likely a fair approximation to ignore the contribution 
from the other images. We will therefore assume that fi°^ (u Q ) / fj, (u ) ^> 1, and thus 
6 d (u ) ~ iold ' s (u ). Consequently, 

0ir P K) = fold ' d K) - fold '>o) (91) 

Since 8 Md ' d (u ) ^ 8 iold,s (u ), we obtain a discontinuous jump in the image centroid as 
the source passes through the fold point u . There is no discontinuous jump if the source 
passes through the cusp (since 6 d mp (0) = 0). Since by equation (25) we have (— qo/2) 1 / 3 = 
-F( — Po/3) 1 ^ 2 , where the minus (respectively, plus) corresponds to uq on the fold arc below 
(respectively, above) the Mi-axis, equations (46) and (47) yield 



<TO = I ~«,,*>t/-F ] - I ~i^y^^ I • (<«> 

where as before the =F signs are for u on the respective fold arcs below and above the tangent 
line to the cusp (i.e., the Mi-axis). Note that bu\ < by equation (40). 

For uq sufficiently close to the cusp at 0, the variable vf{ dominates (w?) 2 . We can then 
approximate the magnitude of the jump as follows: 

1/2 



b 2 u® \ 2 -66m? 
1 + 



c(2ac - b 2 ) ) 2ac - b 2 



-Qbu\ 
2ac - b 2 ' 



(93) 



Therefore, the magnitude of the centroid jump increases as |ii?| 1//2 for \u°\ <C 1, where 
\u*l\ is the horizontal distance of the caustic point u Q from the cusp. Note that no jump 
discontinuity occurs for a source passing through the cusp. 



2.4. Extended Sources Near a Cusp 

2.4-1- Magnification and Image Centroid of Extended Sources 

Consider an extended source with surface brightness S(u). The average surface bright- 
ness is S = {upl)^ 1 j v du S(u), where D is the disc-shaped region of the source and 
= 9^/9 E is the angular source radius 9* in units of E . The magnification of the finite 
source is then 

_ J v duS(u)^u) 
" " f v duS(u) ' (94j 
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Introduce new coordinates by 

it = , (95) 

p* 

where u cn is the center of the source. Note that \u'\ < 1. Then equation (94) simplifies to 

// s = - / du'S N (u')fi(u'), (96) 

J D 

where S N (u) = S(u)/S. In the case of small source sizes 9* <C 9 E , the magnification /x (w) 
of the images not associated with the cusp is a slowly varying function of u over the source, 
which yields 

/4 S ^ M^cn) = /iO,cn- (97) 

The total finite source magnification becomes 

/"t f ot = - / du'SjvKVtotK) = ^ + (98) 



where 



= - f du'^KKcK)- (99) 



For an extended source, the image centroid is given by 

a fs _ J v duS(u)fi tot (u)O d (u) f v du'S N (u')p tot (u')6 cl (u') 

tf i = — — r — r = 7 , I -Will 

J v duS{u)p tot (u) 7Vpf ot 

where the denominator of the middle ratio is TrplS/if ot . For a source with angular radius 
0* <C $Ej the centroid 6q{u) of the images not associated with the cusp varies slowly over 
the source. Consequently, since po(u) also varies slowly over the source, we get 

V0,cn d 0,cn - Po(u cn )0 Q (u cn ) = /i , C n^0,cn, (101) 

where we have defined 0o,cn = Oo(u cn ). Equation (85) then yields 

fs 

Ptot Ptot 

where the contribution from the images affiliated with the cusp is 

fs _ J v duS(u)fi loc (u)O loc (u) _ J v du' S N (u')fi loc (u')O loc (u') 

loc " f v duS(u)p loc (v) *l& ' [ } 

and from the images unrelated to the cusp is 0o,cn- 
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2.4-2. Small Finite Sources on the Axes 



We would like to be able to apply the formulae in the previous section (§2.4.1) to the 
analytic expressions derived for the sources near a cusp (see §2.2.5), to find analytic (or 
semi-analytic) expressions for the finite-source magnification and image centroid of arbitrar- 
ily large sources with arbitrary positions with respect to the cusp. Unfortunately, due to 
the complicated forms for these quantities for arbitrary point source positions, this is effec- 
tively impossible, and numerical methods must be employed. We take this approach in §3. 
However, it is possible to obtain analytic expressions under certain restrictive assumptions. 
Specifically, we now determine the local magnification and local image centroids for relatively 
small uniform and limb-darkened extended sources on the u±- and -u 2 -axes. 

The magnification and image centroid of an extended source are given by equation (98) 
and equation (102), respectively, where the key terms are the local magnification (Eq.[99]) 
and local image centroid (Eq. [103]). The quantitative forms for the magnification and cen- 
troid of a finite sized source depends on the specific form of the surface brightness profile 
S(u). In this section, we derive simplify expressions that involve integrals over arbitrary 
surface brightness profiles. In the next section (§2.4.3), we adopt specific surface bright- 
ness profiles to derive semi-analytic expressions for the centroid and magnification for small 
sources on the fold axes. 

We first consider sources on the Mi-axis, such that V is entirely exterior to the caustic, 
i.e., such that C(u) > for all u in V. Express the center of the source as u cn = (zp*,0), 
where z > 1 with z — 1 corresponding to the boundary of the source meeting the cusp 
and z > 1 yielding a source on the positive ^i-axis away from the cusp. We shall assume 
that for a sufficiently small source /)„ C 1, points u = (ui,u 2 ) inside the source are such 
that «2 ~ and so equation (26) yields 9 2 (u) pa 0. By equations (30) and (31), the local 
magnification at points u inside a sufficiently small uniform source centered at u cn = (zp*, 0) 
can be approximated by p°^(u) = \bui\~ 1 , or, equivalently, 




(104) 



Using equation (104), the finite source local magnification (Eq.[99]) becomes 




1 



Vs N (z), 



(105) 



P*\b\ 



where 




(106) 
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Now, at points u inside the source, equation (28) and equation (38) show that the 
associated local image centroid can be approximated as follows (since 9 2 (u) 0): 



<C'(«) 



^,o) 

C / 



^loc l^cnj H 1, 



(107) 

c / c 

where i = (1,0). Inserting equation (107) into the middle ratio in equation (103) yields a 
simple formula for the local image centroid: 



fs,out 
oc 



where 



\b\c£ 



Qs N (z), 



d-u' 2 S , Ar(-U / L , u' 2 ). 



(108) 



(109) 



For sources on the Mi-axis for which T> is entirely interior to the caustic, i.e., such that 
C(u) < for all u in V, the resulting expressions are quite similar to those for sources 
exterior to the caustic. In particular, the local magnification at points u inside a sufficiently 
small uniform source centered at u cn = (—zp*,0) can be approximated by, 



and therefore 



fs.iri / 



p*\b\ \u[ + z\ ' 
2 



p*\b 



(110) 



111) 



Similarly, the centroid of the three images created when the source is interior to the cusp 
can be shown to be, 



1 + 



2{2ac-b 2 ) 



,0 



Ocn) + P ~ 



1 + 



2(2ac - b 2 ) 



Mil, 



(112) 



which is identical to the centroid for sources exterior to the caustic (Eq.[107]) with the 
exception of the term in square brackets. The finite source image centroid is therefore, 

C n W = 0!ocKn)+ ^ 



H C /4oc m 



1 + 



2(2ac - b 2 ) 



Qs N (z), 



(113) 



Finally, we consider the magnification of small sources on the u 2 axis, for which V is 
entirely exterior to the caustic, i.e., such that C(u) > for all u in V. For sources on the 
■u 2 -axis, the magnification is 



Ur2 
P* 



2/3 



K + P*) 2/r 



(114) 
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The finite source local magnification becomes 

^ ±( * ) = ur) UsAzl (115) 

where 

1 /-i i /VMS) 2 

Ks N (z) = - du[ / dw'^ivK,^). (116) 

We have found similar (semi-) analytic expressions for the centroid of small sources centered 
on the w 2 -axis; however these expressions are generally unwieldy, and thus not very useful 
analytically. Therefore, for the sake of brevity, we will not present them here. 



2.4-3. Uniform and Limb- Darkened Sources on the Axes 

In this section, we will consider two forms for the surface brightness profile. The first, 
and simplest, is a uniform source, where the normalized surface brightness obeys 

S N (U) EE 1. (117) 

Second, we consider the following form for the normalized surface brightness, 



s N (u) = <i-r 



2\ 1/2' 



2 V pI 



(118) 



This form, first introduced by Albrow et al. (1999b), is applicable for microlensing in the 
local group, and is appropriate to a non-uniform stellar source with limb darkening. This 
profile is parameterized by the limb-darkening coefficient T, the value of which is typically 
dependent on wavelength. A key feature of equation (118) is that no net flux is associated 
with the limb darkening term. Note that a uniform source is simply the specific case of 
equation (118) with T = 0. 

Adopting these forms for Sn, we can determine the integral functions Vs N , Qs N , and 
lZs N , which dictate the basic astrometric and photometric behavior of small finite sources 
on the axes of the cusp. The following special case of Vs N (z) (Eq.[106]) is technically handy: 



V n (z) = - T dx^— I 1 X dy(n + 1) [(1 - x 2 ) - y] n , (119) 

where n = k/2 with k > —2 an integer. Using the following identity (e.g., Zwillinger 1996, 
p. 387, Eq. 599), 



v a 2 - y 2 ) n dy = . n \ ,^ 2{n+1/2 \ (120) 
v y ' y 2 (n + 1/2)! ' v ; 
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1 (n+l)\ f 1 , (\-x 2 ) n+1 l 2 

(IX : : . 



equation (119) simplifies to 

V (z) ■ 

ny 1 (" + 1/2)17-1 \x + z_ 

Similarly, we can identify the special case of Qs N (z) (Eq.[109]), 

1 
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dy(n + l)[(l-x 2 )-y] n , 



1 (n + l)\ f 1 , (\-x 2 ) n+l / 2 
dxx- 



xAF(n + l/2)! 



y/(x + zf 



(122) 



and TZ Sn {z) (Eq.[116]), 

1 f 1 

^ J-i 



dxx- 



+ z) 2 / 3 



[X 



1 (n + 1)! 



0F(n + l/2) 



dxx 



[1-x 2 ) 



dy(n+l)[(l-x 2 )-y] n , 



2\n+l/2 



(123) 



For a small uniform source outside the cusp, and on the axis of the cusp (i.e., tangent 
line), equation (105) shows that the local magnification becomes 



us, out/ \ 



p*\b\ 



Mz) 



(124) 



Similarly, from equation (110), the magnification for a small uniform source inside the cusp, 
and on the axis of the cusp, is 



us, in/ \ 
Ploc (*) 



(125) 



By equations (108) and (124), the uniform source local image centroid for source exterior to 
the caustic is 



CrV) = ^oc Kn) + 



P* Qo(z) ^ 
c V (z) h 



(126) 



where 6C(w cn ) = (zp./c,0). For sources interior to the caustic, the corresponding centroid 
is, 



r> us, in/ \ /jin / \ , P* 
O loc ( Z ) = ^locl^cn) + — 



1 + 



2(2ac - b 2 ) 

For a limb darkened source, equation (106) becomes 

Vf N (z) =V (z) + T[V 1/2 (z) -V (z)]. 



Qo(z) , 
Vo(z) 1 - 



(127) 
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Consequently, the finite source local magnification exterior to the cusp is 



Id, out 



( z ) = TTn( p o(z) + r[V 1/2 (z)-V (z)]) 1 
p*\o\ 



and interior to the cusp, 



Id, in/ \ 



P*\b\ 



(V (z) + T[V 1/2 (z)-V (z)}) 



Now, equations (109) and (122) yield that 

Q l i N ( z )-Qo(z) + r[Q 1/2 (z)-Q (z)}. 



(129) 



(130) 
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Thus, using equation (131) we see that the local centroid exterior to the caustic (Eq.[108]) 
becomes 



p„ Qo(z) + r[Q 1/2 (z)-Qo(z)l 



locV cn; ' c Vo(z) + r[V 1/2 (z)-V (z)} 
and interior to the caustic (Eq.[113]), 

b 2 1 Qo(z)+r[Q 1/2 (z)-Q (z)}. 



i + 



2(2ac-b 2 ) 



V (z)+T[V l/2 (z)-V (z)} 



(132) 
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Figure 5 shows the characteristic magnification functions Vo(z) and Vi/ 2 (z), as well as the 
centroid functions Qo(z) and Qi/ 2 (z). Also shown is the fractional difference (Vi/ 2 (z) — 
Vo(z)) /Vo(z) and the absolute difference Qi/ 2 (z) — Qo(z). 

Finally, we can consider the magnification of small sources on the -u 2 -axis. Using equa- 
tion (116), for a uniform source we find that 

n^ N (z) = n (z), 

u r2 \ 2/3 



so 



p* J 



Kq{z). 



US,_L/ \ 

For a limb-darkened source, we obtain, 

K l l N {z) = Ko(z) + T [K 1/2 (z) - K (z)l 

and, thus, 
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u r2 \ 
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(n (z) + T [K 1/2 (z)-K (z)}). 
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Figure 6 shows the functions TZo(z) and 7Zi/ 2 (z), along with their fractional difference 

(n 1/2 (z)-n (z))/n (z). 



-28- 



3. A Worked Example: Binary Lensing Event MACHO-1997-BUL-28 

In this section we numerically calculate the observable microlensing properties for a 
binary-lens cusp caustic crossing. We do this in order to illustrate the photometric and as- 
trometric lensing behavior near a cusp, and to provide an estimate of the magnitudes of the 
effects of finite sources and limb darkening for a general source trajectory. In the next section 
(§4), we verify and explore the accuracy of the analytic formulae for the total magnification 
and centroid derived in the previous sections by comparing them to the numerical results 
obtained here. We calculate the photometric and astrometric behavior for the Galactic bi- 
nary microlensing event MACHO 97-BLG-28, which was densely monitored by the PLANET 
collaboration (Albrow et al. 1999a). The best-fit model for this event has the source crossing 
an isolated cusp of a binary lens; the dense sampling near the cusp crossing resulted in a de- 
termination of not only the dimensionless source size p*, but also limb-darkening coefficients 
in both the I- and V-bands. Combined with an estimate of the angular size 9* of the source 
from its color and apparent magnitude, the measurement of p* yields the angular Einstein 
ring radius, #e = @*/p*- Therefore, the absolute angular scale of the astrometric centroid 
shift is known, and thus, up to an orientation on the sky and subject to small errors in the 
inferred parameters, the astrometric behavior can be essentially completely determined from 
the photometric solution. 



3.1. Formalism and Procedures 

Consider a lens consisting of two point masses located at positions 0ii and i2 , with 
no smoothly distributed matter or external shear. In this case, the dimensionless potential 
(Eq.[4]) is given by, 

i/>(0) = mi In \0 - M | + m 2 In \0 - J)2 |, (136) 

where mi and m 2 are the masses of the two components of the lens in units of the total 
mass. Note that all angles in equation (136) are in units of the total mass of the system. 
Since a = V^, the lens equation (Eq. [3]) becomes, 

— 0ii — 012 

u = °- mi w^et? " m2 l^tr (137) 

Equation (137) is equivalent to a fifth-order polynomial in 0, thus yielding a maximum of 
five images. All of the image positions for a given point on the source plane can be found 
numerically using any standard root finding algorithm. Then the individual magnifications 
can be found using equation (6), and the total magnification and centroid of these images 
are given by equations (7) and (8). 
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To calculate the total magnification and centroid for a finite size source, it is necessary 
to integrate the corresponding point-source quantities over the area of the source (see Eqs. 
94 and 100). This is generally difficult or time-consuming to do directly in the source plane, 
because of the divergent nature of these quantities near the caustic curves; furthermore, it 
requires one to solve equation (137) for each source position. Instead we use the method 
of inverse ray shooting, which is generally easier, quicker, and more robust. We sample the 
image plane uniformly and densely, and use equation (137) to find the source position u{0) 
corresponding to each trial image position 6. The resulting source positions are binned, 
with the magnification of each bin equal to local ratio of the density of rays in to the 
density of rays in u. The result is a map of the magnification //(it), the magnification as 
a function of u. Similarly, one can determine the astrometric deviation by sampling in 0, 
using equation (137) to determine u(0), and then summing at each u the values of 0{u). 
The astrometric deviation at u is then the summed values of 0(u), weighted by the local 
magnification. Thus, one creates two astrometric maps, for each direction. The resolution 
of the resulting maps is given by the size of the bins in the source plane. Since inverse 
ray shooting conserves flux, the maps can be convolved with any source profile to produce 
the finite source photometric and astrometric behavior for arbitrary source size and surface 
brightness profile. 

In this method, the accuracy of a resolution element is limited by Poisson fluctuations 
in the number of rays per bin. The fractional accuracy in the magnification /i k of resolution 
element k is therefore 



where a is the surface of density of rays in the image plane and Q the area of each resolution 
element. For a source size encompassing many resolution elements, the total error is 



where the sum is over all the resolution elements, and Wk is the weight of each resolution 
element. For the simulations of event MACHO 97-BLG-28, we require that the resolution 
of the maps be considerably smaller than the dimensionless source size p*, and we require 
extremely high accuracy in order to explore the subtle effects of limb-darkening. For MACHO 
97-BLG-28, the best-fit model 5 has = 6>*/6> E = 0.0287. We use a resolution element of 



5 We will be focusing on the first "LD1" model of Albrow et al. (1999a), with parameters given in column 
two of their Table 1. This corresponds to the best-fit model assuming a one-parameter limb-darkening law 
and DoPHOT error bars. Although their model LD2 provides a marginally better fit, the differences arc not 
important for this study. 




(138) 




(139) 
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angular size 2 x 10~ 4 6 E ~ 0.007 9* and, thus, fl = 4x 10~ 8 #|. There are approximately 
6.5 x 10 4 resolution elements per source size. We sample the image plane with a density 
of a ~ lO 9 ^ 2 and, hence, (cr M ///)fc = 15% / u fc 1 ^ 2 . Therefore, the total fractional error for 
each source position is always < 0.1%, considerably smaller than any of the effects we are 
considering. 

3.2. Global Astrometric Behavior 

We first analyze the global photometric and astrometric behavior of event MACHO 
97-BLG-28, which is summarized in Figure 7. The best-fit model of Albrow et al. (1999a) 
has the source being lensed by a close-topology binary lens. The topology of a binary-lens is 
specified by the mass ratio, q = 1712/1711 = 0.232 and the projected separation in units of the 
total #e of the system, d = \0i t i — 6i^\ = 0.687. Furthermore, we shall employ a coordinate 
system such that the masses are located on the #i-axis, with the origin equal to the midpoint 
of the two masses and the heavier mass being mi, which is located on the right. Thus, we 
have Oi : i = (d/2, 0) and 6 i2 = {—d/2, 0). For these parameters (d, q), the caustics form three 
separate curves, with the largest (primary) caustic located between the two masses, having 
four cusps, two of which have axes of symmetry that are coincident with the #i-axis. Since 
mi 7^ m 2 , the caustic is not symmetrical about the #2-axis, and there is a relatively isolated 
cusp located approximately at the origin. See Figure 7(a), where we show the caustics and 
source trajectory for the MACHO 97-BLG-28 model. 

For source positions that are large in comparison to d, the lens acts like a single mass 
located at the center of mass of the binary. Thus, for such a close binary lens, two of the 
images behave similarly to those formed by a single lens. This can be seen in Figure 7(b), 
where we plot the critical curves and image trajectories for the event. Similarly, the total 
magnification behaves like the magnification of a single lens well away from closest approach 
to the origin (midpoint of the binary). It is only when the source moves relatively close to 
the lens that the binary nature of the lens is noticeable. See Figure 7(c), where we plot the 
magnification as a function of time for 70 days centered on the maximum magnification of 
the event. Inspection of Figure 7(c) shows that, at the time of the closest approach to the 
origin of the lens, the source passes very close (0.003#e) to the origin. It thus traverses the 
isolated cusp of the primary caustic, resulting in a spike in the magnification (Fig. 7c). This 
also results in large excursions in both components of the image centroid, as can be seen in 
Figure 7(c,d,e), where we plot the two components of 6 C \ as a function of time (panels c and 
d), and also the trajectory of centroid shift 86 c \ (panel e). Note that 6e — 300/ias (Albrow 
et al. 1999a), which sets the scale for 6 C \ and 56 c \. 
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3.3. The Cusp Crossing 

Dense photometric sampling of the event allowed Albrow et al. (1999a) to not only 
determine the global parameters of the lens, but also the source size p* and limb-darkening 
coefficients T in two different photometric bands, Tj = 0.67 (/-band) and T v = 0.87 (V- 
band). See equation (118) for the definition of T. We can therefore explore not only the 
photometric but also the astro-metric behavior of the cusp crossing, including both the finite 
source size and limb darkening. Furthermore, since 9e is approximately known, we can 
determine the absolute scale of these astrometric features and assess their detectability by 
comparing them with the expected accuracy of upcoming interferometers. 

Figure 8 summarizes the astrometric and photometric properties near the cusp cross- 
ing. Figure 8(a) shows the cusp and source trajectory; the cusp point is located at u\ fi = 
—0.009597, and at closest approach, the center of the lens is separated by ~ 0.0066 9e, or 
~ 20% of a source radius, from this point. Figure 8(b) shows the positions of the images as- 
sociated with the cusp at intervals of 15 minutes, along with the parabolic critical curve. The 
light curve near the cusp crossing is shown in Figure 8(c), for a point source, uniform source, 
and limb-darkened sources in I and V. The difference between the point source and uniform 
source is considerable; whereas the difference between the uniform and limb-darkened source 
is smaller (but still significant). In Figure 8(d,e), we show the two components of 6 C \ as 
a function of time for the same source profiles, and in Figure 8(f), we show the trajectory 
of 56 c \. While the nearly instantaneous jumps are smoothed out by the finite source size, 
and also magnitude of the deviations are somewhat suppressed, ~ 200yuas excursions due 
to the cusp are still present. Such variations should be readily detectable using upcoming 
interferometers for sufficiently bright sources. 

We show the path of the centroid 6 C \ near the cusp in detail in Figure 9, for a point 
source, uniform source, and limb darkened sources with Ti = 0.67 and Ty = 0.87. As the 
source crossed the bottom fold, two images appear at the critical curve on the opposite side 
of the axis of symmetry of the cusp from the source, resulting in a large, discontinuous jump 
in 6 C \. Interior to the caustic, there are three images and the centroid moves in the opposite 
direction to the motion of the source, until a pair of images disappear when the source is on 
the top fold, resulting in another discontinuous jump in 6 C \. Again, these images disappear 
on the critical curve on the opposite side of the axis as the source. This motion results in a 
characteristic "swallowtail" trajectory that can also be seen in Figure 3, and is characteristic 
of the centroid of a pointlike source crossing both of the fold arcs associated with a cusp. 
For a sufficiently large source size, as is the case for MACHO 97-BLG-28, this behavior is 
completely washed out, resulting in a smooth trajectory. Not surprisingly, the maximum 
astrometric deviation is also somewhat suppressed. In this case, the maximum difference 
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between the astrometric shifts for point source and finite source amounts to ~ 50/ias. 

In Figure 10, we show the two components of = 0^—6™, the astrometric offset due 
to limb darkening, as a function of time. The are two important points to note. First, the 
deviation is largest just after the source limb first crosses the fold caustic, and just before 
the source limb leaves the caustic entirely. Second, the two components of are not 
correlated, implying the limb-darkening offset is truly two-dimensional. This is opposed to 
the astrometric offset in pure fold caustic crossings, which is essentially one-dimensional, at 
least for simple linear folds (Paper I). This latter point is well-illustrated in Figure 11, where 
we show the trajectory of A0^. The maximum astrometric shift due to limb-darkening 
for MACHO 97-BLG-28 is ~ 35Fyuas. This is similar to the maximum astrometric shift 
due to limb-darkening for the fold caustic-crossing event OGLE 1999-BUL-23 in Paper I 
(~ 60r,uas), once the factor of ~ 2 difference between the value of 9e of the two events 
is taken into account. However in the case of the cusp, the deviation is > 10/ias for a 
considerable fraction of the time the source is resolved, whereas in the case of the fold, the 
deviation is only > 10/ias for the last ~ 10% of the caustic crossing. 



4. Comparison with Analytic Expressions 

In § 2, we derived analytic expressions for the local photometric and astrometric behavior 
of sources near cusp singularities, using the generic form for the cusp mapping, whereas in 
§3, we used numerical techniques to study the behavior of microlensing observables near the 
cusp of a binary-lens. It is important to make contact between the analytic and numerical 
results, and, in particular, to explore the validity and applicability of the analytic results to 
observable lens systems. 

In order to compare the two results, we must specify the coefficients a, b, and c of the 
generic cusp mapping equation (9), that correspond to the cusp that was crossed in the 
binary-lens event MACHO 97-BLG-28. These coefficients are defined in terms of the local 
derivatives of the binary-lens potential at the point on the critical curve corresponding to 
the image of the cusp point. The critical curves of a binary lens can generically be found by 
solving a quartic complex polynomial (Witt 1990). However, for the binary- lens coordinate 
system described in the previous section (masses on the 6*i-axis and origin at the midpoint 
of the lenses), the equation for critical points on the 6\ axis takes the simpler form, 

rm [~\ - c ,i) 2 + ™ 2 (± - 9^ 2 - (± - 9 C ^ 2 (- d - - 6 C ^ 2 = 0. (140) 
For d = 0.687 and q = mi/m 2 = 0.232, the cusp critical point of interest is C = (—0.9488, 0), 
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which is mapped to the cusp caustic point u c = (—0.0096, 0) using the binary-lens equation 
(Eq.[137]). 

The coefficients a, b, and c can be found using the expressions (Eq.[10]) by taking the 
derivatives of the potential (Eq.[136]). We find, for critical points on the ^-axis, that 

mi m i 
(9 cA -d/2) 4 + (^ Cil + rf/2) 4 ' 

(9 cA -d/2) 3 + ^c^ + d^) 3 ' 
l | mi (9 Ctl + d/2) 2 + m 2 (9 C)1 -d/2) 2 
(d/2- 9,, f [dyl ■ 9 rA f 

Inserting equation (140) into equation (143) we find that c = 2 for any critical points on the 
0i-axis. For C = (-0.9488, 0), we find a = 1.694, b = -2.450, c = 2. With these coefficients, 
the local behavior of the cusp is completely specified, and the magnifications and positions 
of the images associated with the cusp can be determined using the expressions derived 
in § 2.2. However, in order to compare with the photometric and astrometric behavior of 
MACHO 97-BLG-28, not only the local behavior of the cusp images must be specified, but 
also the behavior of the images not associated with the cusp. As in § 2.3, we assume that the 
magnification and centroid of the two images unrelated to the cusp can be well-represented 
by a first-order Taylor expansion about the time t c when the source crosses the axis of the 
cusp (i.e., tangent line to the cusp). We find the slope and intercept of these expansions 
that minimizes the differences between the analytic and numerical calculations. Note that 
these parameters (the three slopes and three intercepts) are the only free parameters, the 
photometric and astrometric behavior of the images associated with the cusp are completely 
determined once the coefficients a, b, and c are specified. 

In Figure 12, we show the point-source photometric and astrometric behavior near the 
cusp crossing of MACHO 97-BLG-28, as determined using both the full binary-lens formalism 
and using the analytic expressions. Also shown in Figure 12 are the corresponding predictions 
for a finite source. To calculate the finite source magnification and centroid for the generic 
cusp, we adopt the inverse ray-shooting procedure described in § 3.1, using the cusp mapping 
(Eq. [9]) to predict the source positions. In both the point-source and finite source case, the 
agreement is quite remarkable. We find that the difference in magnification is < 5%, and 
the difference in the centroid is < 5/zas. This is generally smaller than any of the physical 
effects we have considered, including limb-darkening. 
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5. Some Applications 



Generally speaking, the results derived in §2 on the photometric and astrometric mi- 
cro lensing properties near cusps are useful for two basic reasons. First, they are (for the 
most part) analytic, making computations enormously simpler and considerably less time- 
consuming, and making results easier to interpret. This allows for rapid and efficient ex- 
ploration of parameter space; a significant advantage over numerical methods. Second, the 
results are generic; they are applicable to any cusp produced by a gravitational lens. Sim- 
ilarly, they are local, being tied to the global lens in consideration only via the coefficients 
a,b,c of the generic cusp mapping (Eq. [9]). This means that the observables due to the 
cusp crossing can be analyzed in a parametric manner (with the coefficients a, b, c being the 
parameters), without regard to the global (and presumably, non-analytic) properties of the 
lens in question. In this section, we discuss several possible applications of our results that 
take advantage of these properties. 

One issue often encountered in microlensing studies is the question of when the point 
source approximation is valid. This is important, because it is typically easier to calculate 
the observable properties under the assumption of a point-source than when a finite sources 
is considered. For a given source location, the accuracy of the point source approximation 
depends primarily on the distance of the source from the nearest caustic, relative to the size of 
the source, and secondarily on the surface brightness profile of the source. For a source near 
a simple linear fold, there exist semi-analytic expressions for the finite source magnification, 
relative to the point-source magnification, for uniform and limb-darkened sources (see, e.g. 
Paper I). For a source interior to the fold caustic (where there are an additional two images 
due to the fold), the difference in the uniform source magnification from the point source 
magnification is < 5% for separations z > 1.5 , where z is the distance of the center of 
the source from the fold in units of the source radius, and the point-source magnification is 
proportional to 2T 1 / 2 . The effect of limb-darkening is < 1%T in the same range, where T 
is the limb-darkening parameter, which has typical values in the range 0.5 — 0.7 for optical 
wavelengths. Using the results derived in §2.4.3, we can also determine the effect of finite 
sources on the magnification near a cusp. For sources outside the cusp and on the axis of 
the cusp, the point source magnification is fj,°^ = u r \/{p*z) (Eq.[36]), whereas the uniform 
source magnification is /^ out = (u r i/ p*)Vo(z) (Eq.[124]). Thus, the fractional deviation 
from a point source is zVo(z) — 1, independent of and u r \ (at least for small sources). 
Figure 5(a) shows zVo(z) as a function of z. We find that the uniform source magnification 
deviates from the point source magnification by < 5% for z > 2.5. Similarly, one can show 
that the effects of limb-darkening are < l%r for z > 2.5; see equation (129) and Figure 
5(c). For a given distance to the cusp, the effects of finite sources are less severe for sources 
perpendicular to the axis of the cusp than parallel to the axis (see equations [134] and 
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[135]). Thus, one can use the results for sources along the axis of the cusp to determine 
the applicability of the point-source approximation for the magnification given a certain 
photometric accuracy and arbitrary source trajectory relative to the axis of the cusp. It is 
interesting to note that, at a given distance from the caustic, the effects of finite uniform 
and limb-darkened sources are larger for cusps than for folds. 

The remarkable agreement of the analytic predictions for the properties of MACHO 
97-BLG-28 with the exact numerical calculation suggests another application of our results. 
One primary difficulty with modeling Galactic caustic-crossing binary lenses is the fact that 
the magnification is generally non-analytic. Furthermore, in order to calculate the finite 
source magnification, this non-analytic point-source magnification must be integrated over 
the source. This can be quite time consuming, thus making it difficult to explore quickly and 
efficiently the viable region of parameter space. Furthermore, the observables are generally 
not directly related to the canonical parameters. This means that small changes in the pa- 
rameters generally lead to large changes in the observables. Therefore, the parameters must 
be quite densely sampled, further exacerbating the problems with the calculations. Thus, 
finding all viable model fits to well-sampled Galactic binary lensing events has proved quite 
difficult. These problems have been discussed in the context of fold-caustic crossing binary 
lens events in Albrow et al. (1999b), who proposed an elegant and practical solution. Taking 
advantage of the fact that the behavior near a fold caustic is universal, and using analytic 
expressions for the point-source and finite-source magnification, Albrow et al. (1999b) de- 
vised a procedure wherein the behavior near the fold crossing is isolated and fit separately 
from the remainder of the light curve. Again, because the behavior near a fold can be cal- 
culated semi-analytically, fitting this subset of the data is considerably easier and quicker 
than fitting the entire data set. This fit can then be used to constrain the viable regions 
of the parameter space of the global binary-lens fit. We propose that a similar procedure 
might be used to fit cusp-crossing binary-lens events. Consider an observed binary-lens event 
with one well-sampled cusp crossing (e.g. MACHO 97-BLG-28). Extracting the portion of 
the light curve near the cusp crossing, a fit to the generic cusp forms presented in §2 can 
be performed, allowing for images not associated with the cusp. This fit can then be used 
to constrain viable combinations of the parameters a, b, c. One can then search for cusps 
in the space of binary-lens models that satisfy these constraints. Finally, one can search 
for viable fits to the entire data set in this restricted subset of parameter space. Although 
there are undoubtedly nuances in the implementation, it seems likely that this (admittedly 
schematic) procedure should provide a relatively robust and efficient method for fitting many 
cusp-crossing binary-lens events. 

We discuss one final application of our results that may prove useful in the context 
of quasar microlensing. For poorly-sampled microlensed light curves, there exists a rough 
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degeneracy between the effects of the size of the source and the typical angular Einstein 
ring radius (#e) (see, e.g. Wyithe, Webster, & Turner 2000). As discussed by Lewis & Ibata 
(1998), this degeneracy can be broken with astrometric observations, since the scale of the 
motion of the centroid of all the microimages is set by (#e)- Due to the highly stochastic 
nature of the centroid motion for microlensing at high optical depth, in order to use an 
observed data set to determine (6>e), one would have to calculate many realizations of the 
distribution of centroid shifts, and compare these statistically to the observed distribution. 
This procedure is likely to be quite time-consuming. It may be possible to use the results 
presented here to devise a method to obtain a cruder, but considerably less time consuming, 
estimate for (9 E ) given an observed data set. This method makes use of the analytic results 
we obtained for fold caustic crossings. The centroid shift can be quite large when the source 
crosses a fold caustic, when two additional images appear in a position generally unrelated 
to the centroid of the other images. In Paper I, we argued that finite sources only provide a 
small perturbation to the centroid jump when the source crosses a fold caustic, at least for 
sources with angular size significantly smaller than (6>e). For the case of a pure fold caustic 
crossing, it is not possible to provide a prediction for the magnitude of the centroid jump 
from a purely local analysis (Paper I). However, in the case of the source crossing one of 
the two fold arcs abutting a cusp (see Fig.2), the fact that only two images are associated 
with the fold, and one image varies continuously as the source crosses the caustic, implies 
that one can make a definite prediction for the magnitude ly^" 115 ! of the centroid jump for 
a given source trajectory for a generic cusp. In §2.3.2, we found that, for sources crossing a 
fold arc at a horizontal distance from the cusp of the centroid jump is given by, 



where Mj um p = — (2ac — b 2 )/(6b). This expression is valid for \ui\ small, and assuming 
that the total magnification of the images not associated with the cusp is small. Given a 
certain global lens geometry (i.e., shear and convergence), one can obtain an estimate of the 
distribution of the quantity Wj ump by using known techniques to locate the cusps, and taking 
the local derivatives of ip at the corresponding image points to determine the coefficients 
a, b, c. Thus, the distribution of Wj ump can be determined. Similarly, the distribution of u® 
can be estimated once the locations of the cusps are known. Therefore, one can predict 
the distribution of (w^/wjump) 1 ^ 2 for a given lens system. Comparing this to an observed 
distribution of | | , one can constrain (#e)- We note that there are some limitations to 
this method. It is generally only appropriate for sources that are small with respect to (#e)- 
Also, it assumes that the observed centroid shifts are dominated by caustic crossings, and 
that these are well-approximated by the folds abutting cusps. 
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6. Summary and Conclusion 

We have presented a comprehensive, detailed, and quantitative study of gravitational 
lensing near cusp catastrophes, concentrating on the specific regime of microlensing (when 
the individual images are unresolved). We started from a generic polynomial form for the lens 
mapping near a cusp that relates the image positions to the source position. This mapping 
is valid to third order in the image position. The quantitative properties of this mapping 
are determined by the polynomial coefficients, which can be related to local derivatives of 
the projected potential of the lens. Near a cusp, the critical curve is a parabola, which maps 
to the cusped caustic. We find an simple expression for the vertical component of the image 
position 2 , which is a cubic of the form Q\ + p9 2 + q = 0, where p and q are functions of the 
source position u. 

The solutions of the cubic in 62 are characterized by the discriminant D{u) = (p/3) 3 + 
(p/2) 2 . For source positions u outside the caustic, we have D(u) > and, thus, there is 
locally one image. We determined the magnification and location of this image, and showed, 
in particular, that along the axis of the cusp (tangent line to the cusp), the magnification 
scales as u^ 1 , where u is the distance from the cusp point. Perpendicular to the axis, the 
magnification scales as u~ 2 l z . We also determined the image positions and magnifications for 
sources on the caustic, where D{u) = 0. On the caustic, there are two images. One image 
is infinitely magnified and results from the merger/creation of a pair of images, whereas the 
second image has finite magnification, and can be smoothly joined to the single lensed image 
of a source just outside the caustic. For sources inside the caustic (D{u) < 0), we find that 
there are three images. One image has positive parity and diverges as the source approaches 
the top fold caustic abutting the cusp, but can be continuously joined to the non-divergent 
image when the source is on the bottom fold. Similarly, the other positive parity image 
diverges as the source approaches the bottom fold, but can be continuously joined to the 
non-divergent image when the source is on the top fold. The third image has negative parity, 
and diverges as the source approaches either fold caustic. All three images diverge as one 
approaches the cusp point. 

For sources on the axis of the cusp, but interior to the caustic, the total magnification 
of all the images diverges as m^ 1 . We also derived analytic expressions for the centroid of all 
three images created when the source is interior to the caustic. We generalized our results 
beyond the local behavior near the cusp by deriving expressions for the total magnification 
and centroid including the images not associated with the cusp. We further considered 
rectilinear source trajectories, and parameterized this trajectory in order to calculated the 
dependence of the photometric and astrometric behavior on time. In particular, due to the 
presence of the infinitely magnified images when the source is on the fold caustics, but finite 
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magnification outside the caustic, we find that the centroid exhibits a finite, instantaneous 
jump whenever the source crosses one of the two folds abutting the cusp. We present 
a formula for the magnitude of the jump that depends only on the local coefficients of 
the cusp mapping, and the location of the caustic crossing. We note that this magnitude 
decreases monotonically as a function of the horizontal distance between the cusp and where 
the source crosses the cusped caustic curve, so that a source which crosses the cusp point 
exactly exhibits no jump discontinuity. 

Beginning with the appropriate modifications to the formulae for the total magnification 
and centroid for finite sources with arbitrary surface brightness profiles, and combining these 
with the analytic results we obtained for the magnification and centroid for point sources 
near a cusp, we derived semi-analytic expressions for the uniform and limb-darkened finite 
source magnification for small sources on and perpendicular to the axis of the cusp. We 
also derived expressions for the centroid of a small source on the axis for uniform and limb- 
darkened sources. 

In order to illustrate the photometric and astrometric lensing behavior near cusps, and 
to provide order-of-magnitude estimates for the effect of finite sources and limb-darkening 
on these properties, we numerically calculated the total magnification and centroid shift 
for the observed cusp-crossing Galactic binary microlensing event MACHO 97-BLG-28. We 
find that the cusp crossing results in large, 0(9e) centroid shifts, which should be easily 
detectable with upcoming interferometers. We find that limb-darkening induces a deviation 
in the centroid of ~ 35iytas. 

We compared our numerical calculations with the analytic expectation, and found ex- 
cellent agreement. Adjusting only the magnification and centroid of the images unrelated to 
the cusp, and adopting the coefficients appropriate to the cusp of MACHO 97-BLG-28, we 
find that our analytic formulae predict the magnification to < 5% and the centroid to 5yuas, 
for positions within ~ 2 source radii of the cusp. 

Finally, we suggested several applications of our results to both Galactic and cosmo- 
logical microlensing applications. We suggest that one can use our analytic expressions to 
determine the applicability of the point-source approximation for sources near a cusp. In 
particular, we note that the finite source magnification deviates from the point-source mag- 
nification by < 5% for sources separated by > 2.5 source radii. We also discussed how 
the local and generic behavior of the cusp can be used to simplify the fitting procedure for 
cusp-crossing events. Lastly, we outlined a method by which the typical angular Einstein 
ring radius of the perturbing microlenses of a macrolensed quasar might be estimated using 
measurements of the jump in the centroid that occurs when the source crosses a fold, making 
use of the analytic expression for the magnitude of the jump derived here. 
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Despite their apparent diversity, the mathematical underpinning of all gravitational 
lenses is identical. In particular, all lenses exhibit only two types of stable singularities: 
folds and cusps. In Paper I, we studied gravitational microlensing near fold caustics; here 
we have focussed on microlensing near cusp caustics. A generic form of the mapping from 
source to image plane near these each of these types of caustics can be found, and used to 
derive mostly analytic expressions for the photometric and astrometric behavior near folds 
and cusps. These expressions can be used to predict the behavior near all stable caustics 
of gravitational lenses, and applied to a diverse set of microlensing phenomena, including 
Galactic binary lenses and cosmological microlensing. 

We are especially thankful to Meredith Houlton for meticulously checking our calcu- 
lations. B.S.G. was supported in part by NASA through a Hubble Fellowship grant from 
the Space Telescope Science Institute, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under NASA contract NAS5-26555. A. P. was supported in 
part by an Alfred P. Sloan Research Fellowship and NSF Career grant DMS-98-96274. 
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Fig. 1. — Illustration of the basic properties of photometric and astrometric microlensing 
near a cusp, for a source trajectory which passes near, but exterior to, the cusp point, (a) 
The heavy solid line is the caustic, the shaded area indicates the region interior to the caustic, 
where three images are formed; one image occurs for sources outside the bounded shaded 
region. The dotted line is the source trajectory. Each sloped line connects a source position 
u to the point u* on the caustic. The Mi-difference between u* and u is the horizontal 
distance A-u^ = \ui — u* cl \. The magnification of an image is oc (Au^) -1 . See text, (b) 
The solid line is the parabolic critical curve, whereas the dotted line denotes the image's 
trajectory, (c) The magnification of the image is shown in (b) as a function of time. The 
solid line is the total magnification, (d) The ci; i-component of the centroid of all the images 
as a function of time, (e) The fl^-component of the centroid as a function of time, (f) The 
solid line is the critical curve, whereas the dotted line shows the path of the centroid of all 
the images. 
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Fig. 2. — Local gravitational lensing for source passing near a cusp; same as Figure 1, except 
for a trajectory which crosses one fold caustic. The dotted lines in (b) show the image 
path for the portion of the source trajectory in (a) outside the caustic, while the dashed 
and dashed-dot lines are the paths for the other two local images. Panel (c) depicts the 
magnifications of the three local images in (b); the solid line is the total magnification of the 
three images. The dotted path in (f) is the centroid trajectory for the images. 
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Fig. 3. — Same as Figure 1 and 2, except for a trajectory which crosses both upper and 
lower fold caustics. 
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Fig. 4. — Contour plots of the magnification of the various images as a function of source 
position, u. In each panel, contours are in the interval log(/i) = —1, 1, in logarithmic steps 
of 0.1, and the dashed line is the caustic, (a) Contours of magnification of the non-divergent, 
positive-parity image, as a function of u. Outside the caustic, this corresponds to fi°^(u). 
Inside the caustic, this corresponds to ^(u) for u 2 > 0, and n 2 n (u) for u 2 < 0. (b) Contours 
of the divergent, positive-parity image, which corresponds to fJ^(u) for u 2 < 0, and fi 2 n ( u ) f° r 
u 2 > 0. (c) Contours of the divergent, negative-parity image, which corresponds to fj,f(u). 
(d) Contours of the total magnification. See text for the definitions of fi^, [i™, /i 2 n , and 
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Fig. 5. — Basic functions which describe the photometric and astrometric properties of finite 
sources on the i^-axis (the axis of symmetry), as a function of the distance z in units of 
the dimensionless source size p*. (a) The basic functions for the photometric behavior. The 
solid line shows zVo(z), whereas the dotted line shows zVx^iz). The fractional deviation of 
a uniform source from the point source magnification is zVo(z) — 1. (b) The basic functions 
for the astrometric behavior. The solid line shows Qo(z); the dotted line shows Qi/2(z). The 
astrometric deviation from a point source is proportional to <2o for a uniform source, (c) The 
basic functions describing the photometric ((Vi/2(z) — Vq{z)) /Vq{z)\ solid) and astrometric 
(Qi/2(z) — Qo(z); dotted) behavior of limb-darkened sources. 
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Fig. 6. — Basic functions which describe the photometric properties of finite sources on the 
■u 2 -axis (perpendicular to the axis of symmetry of the cusp), as a function of the distance z 
in units of the dimensionless source size p*. (a) The solid line shows z 2 ^TZq(z), whereas the 
dotted line shows z 2 ^TZi/2{z). The fractional deviation of a uniform source from the point 
source magnification is z 2 ^TZ (z) — 1. (b) The basic function describing the photometric 
behavior of limb-darkened sources, (1Z 1 / 2 (z) — 1Z (z))/1Z (z). 
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Fig. 7. — The global photometric and astrometric behavior of the best-fit model for the cusp- 
crossing Galactic binary-lens microlensing event MACHO 97-BLG-28. (a) The shaded areas 
indicate the regions interior to the binary-lens caustic, the solid line shows the trajectory 
of the center of the source, whereas the dotted lines show the width of the source. The x's 
denote the positions of the two masses of the lens; the heavier mass is on the right, (b) The 
thin, dotted lines show the critical curves. The dots show the positions of the images at 
fixed intervals of 1 day. (c) The light curve (magnification as a function of time) for a point 
source (solid line) and uniform source (dotted line). (d,e) The two components of cl as a 
function of time. The solid line is for a point source, dotted line for a uniform source, (f) 
The centroid 50 d of all the images relative to the position of the source. 
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Fig. 8. — The photometric and astrometric behavior of MACHO 97-BLG-28 near the cusp 
crossing, (a) The shaded area is shows the region interior to the caustic. The solid line shows 
the source trajectory, and the dotted line shows the width of the source, (b) The dotted line 
shows the critical curve, the circles are the positions of the images at fixed intervals of 15 
minutes. The size of each circle is proportional to the logarithm of the magnification of the 
image, (c) The magnification near the cusp crossing as a function of time. The solid line is 
for a point source, dotted line for a uniform source, the short-dashed line for limb-darkened 
source in the /-band, and long-dashed line for the \^-band. (d,e) The two components of 
the centroid C \ as a function of time. Line types are the same as (c). (f) The path of the 
centroid shift 50 c \. 
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Fig. 9. — Detail of the path of the centroid 6 C \ of MACHO 97-BLG-28 during a three-day span 
centered on the cusp crossing. The solid line is for a point source, dotted line for a uniform 
source, the short-dashed line for limb-darkened source in the /-band, and long-dashed line 
for the V^-band. 
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Fig. 10. — (a,b) The two components of the astrometric offset due to limb-darkening A6 cl = 
O 1 ^ — 0"f relative to a uniform source as a function of time. The short-dashed line is for 
the /-band, whereas the long-dashed line is for the \^-band. The two components are (a) 
parallel and (b) perpendicular to the axis of the cusp (i.e., tangent line to the cusp), which 
is coincident with the binary axis. 
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Fig. 11. — The path of the astrometric offset due to limb-darkening A6 cl = 6 d — 0™ relative 
to a uniform source. 



3. 



Q5 



-50 
-100 
-150 
-200 

-250 



-300 



1 1 1 1 1 1 1 


n 1 1 1 l ^-^ss»? 


r (c) , 
i i i i i i i 


r 7 i 

i i i i i i i 



894 



895 896 
HJD-2449719 



897 



Fig. 12. — A comparison between the photometric and astro metric properties of the event 
MACHO 97-BLG-28 near the cusp crossing computed using the full binary-lens formalism, 
and using the generic cusp mapping. The coefficients of the cusp mapping were determined 
from the local derivatives of the potential at the position of the image of the binary-lens 
cusp. See text, (a) The light solid line shows the point-source magnification as a function 
of time calculated using the binary-lens formalism, the light dotted line (barely visible) is 
using the analytic expressions derived here. The heavy solid line shows the uniform-source 
magnification calculated using the full binary-lens equation, the short-dashed line is using 
the generic cusp mapping. For comparison, the long-dashed line is the magnification for the 
/-band. (b,c) The two components of the centroid d . Line types are the same as in panel 
(a). 



